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ABSTRACT 



Distributions of kinetic energy by latitude and wave number were 
obtained from harmonic analyses of hemispherical weather charts. 
Correlations between daily changes in these kinetic energies were 
also computed. 

The writers wish to express their gratitude to Professor George 
J* Haltiner, U. S. Naval Postgraduate School, and Commander Paul 
Wolff, Officer-in-charge of NANWEP , for their suggestion of the topic 
and continued assistance and encouragement. 



ii 



TABLE OF CONTENTS 



Section 



Title 



1 * Introduction 

2 Development of kinetic energy 

equations 

3 Correlation coefficients 

4 Kinetic energy distribution 

5 Correlation of kinetic energy 

changes 

6 Conclusions 

7 Recommendations 

8 Illustrations 



Page 

1 

2 

5 

6 

11 

16 

18 

19 



iii 



Fig- 

1 

2 

3 

4 

5 

6 

7 

3 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 



Page 

19 

19 

20 

20 

21 

21 

22 

22 

23 

23 

24 

24 

25 

25 

26 

26 

26 

27 

27 

27 

27 

23 

28 

28 



LIST OF ILLUSTRATIONS 



Title 



Total KE (U) , total KE(u), total KE(v), and 
total E 



KE(u) 
KE(u) 
KE(u) 
KE(u) 
KE(u) 
KE(u) 
KE(u) 
KE(u) 
KE(u) 
KE(u) 
KE(u) 
E and 
KE(U) 
KE (u) 
KE(v) 
KE(v) 
KE (u) 
KE(u) 
KE (u) 
KE(u) 
KE (u) 
KE(u) 
KE (u) 



and 


KE(v) 


VS 


n 


for 


17° 


. 5N. 


and 


KE (v) 


vs 


n 


for 


22° 


. 5N. 


and 


KE(v) 


vs 


n 


for 


27° 


. 5N. 


and 


KE (v) 


vs 


n 


for 


32° 


. 5N. 


and 


KE (v) 


vs 


n 


for 


37° 


. 5N. 


and 


KE(v) 


vs 


n 


for 


42° 


. 5N. 


and 


KE (v) 


vs 


n 


for 


47° 


. 5N. 


and 


KE (v) 


vs 


n 


for 


52° 


. 5N. 


and 


KE (v) 


vs 


n 


for 


57° 


. 5N. 


and 


KE (v) 


vs 


n 


for 


62° 


. 5N. 


and 


KE (v) 


vs 


n 


for 


67° 


. 5N. 



E per unit area vs 0 

and KE (u) / KE(v) all n vs 0 

all n vs 0 

all n vs 0 

all n vs 0, different scale 

and KE (v) n : 1 vs 0 

and KE(v) n = 2 vs 0 

and KE(v) n = 3 vs 0 

and KE(v) n = 4 vs 0 

and KE(v) n : 5 vs 0 

and KE(v) n ■ 6 vs 0 

and KE(v) n = 7 vs 0 



iv 



25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 



Page 

28 

29 

29 

29 

29 

30 

30 

30 

30 

31 

31 

31 

31 

32 

32 

32 

32 



Title 



KE (u) and KE(v) n 
KE(u) n = 9 vs 0 
KE(u) n = 10 vs 0 
KE(u) n - 11 vs 0 
KE(u) n = 12 vs 0 
KE(v) n 5 1 vs 0 
KE(v) n = 2 vs 0 
KE(v) n * 3 vs 0 
KE(v) n = 4 vs 0 
KE(v) n 8 5 vs 0 
KE(v) n : 6 vs 0 
KE(v) n - 7 vs 0 
KE(v) n = 8 vs 0 
KE(v) n - 9 vs 0 
KE(v) n 8 10 vs 0 
KE(v) n : 11 vs 0 
KE(v) n “ 12 vs 0 



8 vs 0 



v 



1. Introduction 



In the Monthly Weather Review of July 1958, Commander Paul Wolff, 

U. S. Navy, prescribed a method for improving numerical weather fore- 
casts. His development decreased the errors due to faulty prediction 
of the movements of very long wave components of the pressure patterns. 

At the conclusion of his article. Commander Wolff suggested that 
empirical methods currently employed could be replaced by terms in 
equations describing the mechanism by which energy is transferred to 
and from long waves. 

The purposes of this paper are to show the distribution of kinetic 
energy by latitude and wave number and to correlate daily changes in 
the kinetic energy of various wave numbers and latitudes with one 
another. 

For this research, data for 120 days during the winter of 1947-48 
were used. These data, furnished by the NAVAL WEATHER RESEARCH FACILITY, 
consisted of mean 500-mb heights for each latitude at 5° intervals 
from 15° N. to 70° N. , as well as amplitudes and phase angles for 
wave numbers 1 through 18 for these latitudes. Because of time re- 
strictions, wave numbers 13 through 18 were neglected in our computa- 
tions. 
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2. Development of Kinetic Energy Equations 

Fourier series were employed to represent the height of the 
500-mb surface in order to investigate the distribution of kinetic 
energy. The form used for each latitude was 

z *= z / z 

whe re 

z'= J: cos 

z = height of 500-mb surface 

z : mean height of 500-mb surface 

A = amplitude of wave number n 

— phase angle in radians with respect to Greenwich 
meridian 

0 - longitude in radians 

The quantities z, z, A., , and J n are all functions of latitude. 

Since wind data are relatively scarce, especially over oceans, 
the geostrophic approximation was used to obtain the wind velocity. 
In accordance with this approximation, the wind components may be 
written as 







(2) 


u = 


- JL J_ i! 

■f A J Q 




u = 


- A ± 




v = 

where 


3.1 Jb’ 

6L cos (p r ) Q- 




U / u 


- total zonal 


component of the wind velocity 


V 


meridional 


component of the wind velocity 




- latitude in 


radians 



g : gravity 

a : radius of earth 

The kinetic energies per unit mass of the mean zonal flow, the 
zonal perturbation wind component, and the meridional wind component 
are 1/2 U' , 1/2 u~ , and 1/2 v 2 . 

Utilizing equations (1) and (2), we obtain 




(3) 



0-> 





1/2 v 4 



- J- / 3 V_1 . 

2- \ Z ST. tin (f / ,v c 





i. 



Multiplying by (' ^ ^ and integrating equations (3) from 

0 to 2TT , we obtain kinetic energies per unit cross-section area for 
the latitude circle. Since the variation of the density is negligible 
compared to velocity, a mean value of density will be assumed. 

The resulting kinetic energies expressed in terms of finite differences 
are as follows: , 



n ' V _L l ilc Si 

U~nJ «- fr/li,)'' 




Subscripts 1 and 2 refer to two adjacent latitude bands. 

V = 1/2 (<P,+ <Pz) 

z - mean height 

^/36 = latitudinal increment 

A common factor has been omitted and values quoted subsequently 



are relative kinetic energies. To obtain values of the kinetic 



energies in cgs units, it is merely necessary to multiply the relative 
kinetic energies by 0.145 x 10^ The following quantities have 
been evaluated on the CDC Model 1604 computer for the relative kinetic 
energies of the mean zonal flow and the perturbation components for 
each wave number: 

KE(U) = k(t) (S x - I.) 

KE( U ) = k «Y— (k - k) }/%k'~ 

KE(v) ■ k 

where 

k (fi) z cos (p /sin‘<^> 
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3. 



Correlation Coefficients 



Correlation coefficients were computed between the daily 
changes in relative kinetic energy of each component originally 
calculated, as well as the sums of these components for each wave 
number and each latitude band. Again the Control Data Corporation 
Model 1604 computer was utilized because of its high speed and large 
memory capacity. 

The product-moment correlation coefficient, r , for relating 
two variables x and y was employed. 

r = covariance (x,y) /g~ r 

* <3 



whe re 



covariance (x,y) s 7 ^ X ; yc — Xj 
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4. Kinetic energy distribution 



In order to appreciate the energy distribution presented below, 
it must be remembered that the kinetic energy was computed using the 
geostrophic approximation of the real wind, which is not entirely 
valid, especially at very low latitudes. The latitude range consider- 
ed is from 15°U .to 70°N. The daily kinetic energies for the 120 days 
were computed for eleven latitude bands, each 5° in width. The units 
used here are arbitrary and do not represent absolute values, but 
illustrate the relative magnitudes of the various components of 
kinetic energy as they are distributed geographically. In order to 
yield absolute kinetic energy in cgs units, the relative magnitudes 
must be multiplied by a factor of 0.145 x 10^ . Occasionally, the 

scale on certain graphs is changed to emphasize salient features. A 
note to this effect is addended so that an erroneous conception of 
relative magnitudes will be avoided. 

The following notation was used in the preparation of figures. 

E ■ total kinetic energy 

KE(U) = mean zonal kinetic energy 

KE(u) s zonal component of the perturbation kinetic energy 
for each wave number 

KE(v) = meridional component of the perturbation kinetic 
energy for each wave number 
n = wave number 

/ = latitude 

The wave numbers, n, range from 1 through 12, and the latitude 
bands discussed are represented by mean latitudes of the bands as 



follows : 
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h = 


17°. 5N 




22°. 5N 


/':= 27°. 5N 


1 1 


32°. 5N 


(fir 


37°. 5N 


42°. 5N 


*7 = 


47°. 5N 


& - 


52°. 5N 


ft, - 57°. 5N 


ft- r 


62°. 5N 


<a,= 


67°. 5N 





Figure 1 illustrates the hemispheric distribution of total kinetic 
energy at 500*mb over the latitudes of 15°N. through 70°N. The total 
kinetic energy, E, is composed of total mean zonal kinetic energy, 
KE(U), total zonal component of the perturbation kinetic energy, 

KE(u), and total meridional component of the perturbation kinetic 
energy, KE(v). 

It is seen from this figure that the mean zonal contribution is 
slightly less than the KE(u) contribution. Of the perturbation com- 
ponents, the contribution of the zonal perturbation, K2(u), is about 
three times as great as the meridional perturbation, KE(v). 

Figures 2-12 depict the relative magnitudes of KE(u) and KE(v) 
contributions of wave numbers 1 through 12. Each figure represents 
the spectrum at a given latitude band. It should be noted that 
these spectra are discrete and have meaning only for integral values 
of wave number. The lines connecting points have been drawn to aid 
visual tracing of fluctuations. Also note that the mean zonal con- 
tribution, KE(U), is not dependent on wave number. 

In summary of figures 2-12, we note several significant features. 
Notably, the spectra reveal that at very low latitudes the maximum 
kinetic energy is contained in wave number 2 of the KE(u) component, 
corresponding to a very long wave length of 180°longitude. The 
maximum in the KE(u) component shifts toward higher wave number as 
latitude is increased, until, at 42P. 5N, wave number 4 is a maximum 
contributor. This indicates that the largest single contribution of 



kinetic energy at this particular middle latitude is associated with 
a wave length of 90°. The maximum shifts back to wave number 2 as 
latitude is further increased, although wave number 4 continues to 
be significant. Above wave number 7 the KE(u) component is relatively 
small at all latitudes. 

The KE(v) contribution is comparatively small at lower latitudes, 
but begins to become significant above 42°. 5N. At these higher 
latitudes, the KE(v) contribution is insignificant in wave number 1, 
but compares in magnitude to the KE(u) in higher wave numbers and 
contributes more to the total energy. 

The profile of the total kinetic energy, E, is as depicted in 
figure 13. The initially high value of kinetic energy at $ 
is undoubtedly a reflection of the trades. The minimum at 0*. cor- 
responds to the winter position of the subtropical high pressure 
belt near 22°. 5N. The maximum at f/W shows the latitudinal position 
of the westerly jet to be around 37°. 5N. Again it must be noted that 
the values at the various latitude bands are discrete. The values 
of total kinetic energy decrease steadily as latitude increases 
northward of 37°. 5N. This is partially a result of the decrease 
of the area of a latitudinal band approaching the pole. The dotted 
line in figure 13 shows the total kinetic energy per unit area. 

It is apparent that the kinetic energy decrease at high latitudes is 
therefore primarily a reflection of a decreased wind intensity. 

The contribution of KE(U) is small at latitudes above , but 
contributes significantly to the total kinetic energy at lower lati- 
tudes except at the subtropical high pressure belt, . (See 

figure 14). 
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The profile of the total KE(u), figure 15, illustrates an overall 
decrease with increasing latitude above , whereas the total KE(v) 
given in figure 16, increases up to 0 / and then decreases. The KE(v) 
profile is also shown with the magnitude scale expanded in figure 17 
to emphasize the latitudinal distribution of KE(v). 

The spectra of pertubation kinetic energy are shown in figures 
2-12. It should be noted that the maximum energy is located at 0t 
for n 2 2, as seen in figure 2, followed by a secondary maximum at (#>- 
for n = 4, as seen in figure 7. The ‘profiles are fairly uniform 
at n 2 4 and above, indicating a consistent distribution for the higher 
wave numbers. The very long waves, n 2 1, 2, 3, have higher values 
at low latitudes and a somewhat erratic distribution. One conclusion 
that may be drawn from this set of graphs is that the long waves, 1-3, 
contribute most of the kinetic energy at < 0 ^ . Figure 14 shows 

that this contribution is largely due to the perturbation kinetic _ 
energy. 

Figures 18-41 show the spectra of perturbation kinetic energy 
contributions. Note that an enlarged scale is used beginning with 
figure 30 for the KE(v) profiles. 

In discussing the total zonal perturbation and mean zonal kinetic 
energy contributions, it was pointed out that these components 
comprised the greatest part of the total kinetic energy in the longer 
waves (n = 1-6) at the lower latitudes ( and below). The KE(v) 
contribution is not significant below 0g . It was also pointed out 
that perturbation kinetic energy is not particularly significant 
at any of the latitudes under consideration in short waves (n = 7-12). 
Therefore, in discussing the perturbation profiles we shall confine 
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ourselves to low wave numbers. Figure 18 shows that the KE(v) 
contribution is negligible for n = 1. Figure 30 shows this component 
with an enlarged magnitude scale. Although the value of KE(v) 
for n = 1 increases with latitude, it is two orders of magnitude 
smaller than the corresponding KE(u) . Figures 19-25 show that the 
KE(v) contribution approaches and sometimes equals the KE(u) con- 
tribution at high latitudes for wave numbers 2-8* 

The KE(v) components are expanded in figures 30-41, as a matter 
of interest, to show the latitudinal distribution of this relatively 
small contributor to the total kinetic energy. 
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5. Correlation of Kinetic Energy Changes 

The correlation of 24-hour kinetic energy changes over the 120- 
day data period was made complete as practicable with the aid of the 
CDC Model 1604 computer. The following items were correlated with 
one another: 

1. Changes in total kinetic energy (all wave numbers, all lati- 
tudes) , 

2. Changes in total kinetic energy at each latitude for all wave 
numbers. 

3. Changes in total kinetic energy for each wave number over all 
latitudes. 

4. Changes in total kinetic energy at each latitude for each 
wave number. 

5. Changes in total irean zonal kinetic energy (all latitudes). 

6. Changes in total mean zonal kinetic energy at each latitude. 

7. Changes in total zonal perturbation kinetic energy (all 

wave numbers, all latitudes). 

8. Changes in zonal perturbation kinetic energy at each latitude 
for all wave numbers. 

9. Changes in zonal perturbation kinetic energy in each wave 
number over all latitudes. 

10. Changes in total meridional perturbation kinetic energy 
(all wave numbers, all latitudes). 

11. Changes in meridional perturbation kinetic energy at each 
latitude for all wave numbers. 

12. Changes in meridional perturbation kinetic energy in each 
wave number over all latitudes. 
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The program was designed to extract only correlation coefficients 
equal to or greater than 0.4. A large number of high positive 
correlation coefficients was obtained, but time limitations have pre- 
vented a complete analysis of the results. The complete printout 
of the correlation coefficients with instructions and identifier table 
for locating the items correlated will be left with Professor George 
Haltiner, U. S. Naval Postgraduate School, Monterey, California. 

Twenty- two negative correlation coefficients in the range, -0.40 
to -0.47, were obtained and these are presented in Table I below. 

It was felt that a negative correlation of kinetic energy changes would 
better indicate the existence of kinetic energy transfer relationships. 
The notation is the same as that used in Section 4 with the addition 
of the symbols JE. > which represent the sum of the indicated 

kinetic energy component over all latitudes and all wave numbers, 
respectively. 

Inspection of these correlations shows the existence of change 
in the total meridional perturbation component at latitude 4T 3 . 5N. 
occurring simultaneously with an opposite change in the (1) mean zonal 
component, (2) total zonal perturbation component, and (3) total 
kinetic energy at latitude 22°. 5N. This negative correlation between 
the total meridional component at 47?5N. and the total kinetic energy 
at 22°. 5N. is a result of the negative correlations between the former 
and total kinetic energy changes in all wave numbers except n “ 2 
and n = 5. 

It is of interest that latitude 22°. 5N. approximately corresponds 
to the position of the subtropical high pressure belt during the winter 
months, whereas latitude 47°. 5N. corresponds to the position of the 
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polar front. At the 500-mb level, which approximates the level of 
non-divergence, the average vertical air mass transport at 47°. 5N. 
would be upward, whereas at 22?5N. it would be downward. 
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TABLE I 



£ KE(v) 


at 


470.5 


VS 


KE(u) 


at 17°. 5 for 


n = 3 




£kE(v) 


at 


47°. 5 


vs 


KE(U) 


at 22°. 5 






£ KE (v) 


at 


47°. 5 


vs 


KE (v) 


at 22°. 5 for 


n = 3 




£kE(v) 


at 

z. 


47°. 5 


vs 


E at 


22°. 5 for n = 


1 




£ KE (v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n = 


3 




£. KE(v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n = 


4 




£ KE(v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n = 


6 




£ KE(v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n = 


7 




£ KE (v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n = 


8 




£ KE(v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n = 


9 




£ KE(v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n = 


10 




£ KE (v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n « 
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£ KE(v) 


at 


47°. 5 


vs 


E at 


22°. 5 for n = 


12 




£ KE(v) 


at 


47°. 5 


vs 


Z KE(u) at 22°. 5 






£ KE(v) 


at 


47°. 5 


vs 


£e at 22°. 5 






£ KE(v) 


at 


47°. 5 


vs 


KE(u) 


at 57?5 for n 


= 1 




KE(v) 


at 


32°. 5 


for 


n = 4 


vs KE(u) at 


32°. 5 


for n 5 10 


KE(v) 


at 


32°. 5 


for 


n = 4 


vs KE (u) at 


27°. 5 


for n * 10 


KE(v) 


at 


32°. 5 


for 


n = 4 


<P 

vs Z KE(u) for n = 


10 


KE(v) 


at 


32°. 5 


for 


n 8 4 


vs KE (u) at 


32°. 55 


for n : 7 


KE (u) 


at 


27°. 5 


for 


n = 4 


vs KE(u) at 37°. 5 for n = 3 


KE(u) 


at 


22°. 5 


for 


n - 6 


vs E at 47°. 5 


for n 
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Between these latitudes there is a net air mass transport 
southward above 500 mb and a net transport northward below 500 mb. 
With the meridional circulation pattern in mind, it is reasonable to 
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assume that there is a kinetic energy interchange between these 
latitude belts. Figure 17 shows that the maximum in the latitudinal 
distribution of the total meridional kinetic energy component exists 
at latitude 4T 3 . 5 N. It is our conjecture that changes in the total 
meridional perturbation kinetic energy component in the region of the 
polar front are caused by ; or result in ; opposite changes in the total 
kinetic energy in the belt of the subtropical high. 

The negative correlation between the changes in the meridional 
perturbation component in n = 4 at latitude 3?. 5 N. and the changes 
in the zonal perturbation component in n 5 7 at the same latitude 
suggest a possible energy exchange between these different wave 
numbers. The negative correlation between the changes in zonal 
perturbation components at latitude 27^5 N. in n = 4 and latitude 
3f.5 N. for n = 3 indicates not only an interchange between wave 
numbers but a possible latitudinal interchange. The negative 
correlations between components in the higher wave numbers would 
not contribute appreciably to the total kinetic energy changes at 
530 mb. 
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6 . 



Conclusions 



In summary of the kinetic energy distributions, it is seen that 
total kinetic energy is high in the trades region, low in the belt of 
subtropical high pressure, a maximum at mid-latitudes near the westerly 
jet, and thereafter decreasing to a minimum approaching the pole, 

a; 

The mean zonal kinetic energy contribution is greatest at mid-latitudes 
with a secondary maximum at the trades, and minima in the doldrums and 
high latitudes. Approximately half of the total hemispheric kinetic 
energy at 500 mb is associated with the mean zonal wind. 

The rest of the total kinetic energy is perturbation kinetic 
energy, which resides for the most part in the longer waves, n z 1-7. 
Both the zonal and meridional components vary with wave number as 
shown earlier. At low latitudes only the zonal perturbation component 
is significant. At high latitudes, the meridional perturbation 
component becomes almost as important as the zonal, except for wave 
number 1. With these points in mind, treating both components of 
perturbation kinetic energy together, we see that wave number 2 con- 
tains the maximum perturbation kinetic energy at low latitudes. 

As latitude is increased the maximum shifts from n = 2 to n = 4 at 
mid-latitude, 47^5 N. Above this latitude the perturbation contribu- 
tion at wave number 2 becomes equal to and then greater than that 
at n s 4. This is not to imply that these wave numbers are the only 
significant ones, they are merely the largest. The combined contri- 
butions to the total kinetic energy from the other wave numbers far 
outweigh the effect of these single maxima. 

On the basis of the few negative correlations considered, it is 
seen that exchanges of kinetic energy exist between different wave 



numbers at the same latitude and between various wave numbers at 
different latitudes. A more detailed quantitative study of these inter- 
changes should prove very useful. 



7 . Recommendations 



A complete analysis of the positive correlation coefficients 
computed for the daily changes in kinetic energy is recommended as 
future investigation. 
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